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^ ■ Abstract 



We explore the first stage of the Symanzik improvement program for lattice 
Dirac fermions, namely the construction of doubler-free, highly improved 
' classical actions on isotropic as well as anisotropic lattices (where the 

temporal lattice spacing, at, is smaller than the spatial one). Using field 
transformations to eliminate doublers, we derive the previously presented 
isotropic D234 action with 0{a?) errors, as well as anisotropic D234 
actions with 0{a'^) or 0{a^,a^) errors. Besides allowing the simulation of 
heavy quarks within a relativistic framework, anisotropic lattices alleviate 
potential problems due to unphysical branches of the quark dispersion 
relation (which are generic to improved actions), facilitate studies of lattice 
thermodynamics, and allow accurate mass determinations for particles 
with bad signal/noise properties, like glueballs and P-state mesons. We 
also show how field transformations can be used to completely eliminate 
unphysical branches of the dispersion relation. Finally, we briefly discuss 
future steps in the improvement program. 



1 Introduction 



Lattice calculations suffer from scaling errors, or lattice artifacts, that typically decrease 
like some power, aP, of the lattice spacing (ignoring logarithmic corrections). Continuum 
results are obtained as a — >^ 0, but the cost of a realistic simulation of QCD, for example, 
grows like like some large power, a~^, of the inverse lattice spacing (a; is at least 6, 
but could even be 10 or more [1]). It is therefore extremely important to find highly 
improved actions; they will give accurate results on much coarser lattices. 

Classical field theory estimates suggest that eliminating errors through order and 
maybe o;^ allows one to model the properties of a smooth bump with errors of a few 
percent to a fraction of a percent by using a lattice with 3-6 grid points per diameter 
of the object in each direction. For hadrons this means that spatial lattices of spacing 
a = 0.2 — 0.4 fm might suffice for improved actions (whereas a = 0.05 — 0.1 fm are 
typically used for unimproved actions). Even considering the computational overhead 
due to the more complicated form of improved actions, it is clear that being able to work 
on coarse lattices would save many orders of magnitude in CPU time. 

For pure glue it has already been demonstrated [2] that this is possible. In this paper 
we describe some of the steps that are necessary to extend these savings to the more 
difficult problem of lattice quark actions. Our considerations are mainly classical, but 
we will outline where quantum effects seem to play a role and how to take them into 
account. The approach we would like to follow, pioneered by Symanzik [3, 4, 5, 6, 7], is 
to try to improve lattice actions and fields to some finite order in a, like or a^. 

For asymptotically free theories, such as QCD, the terms in the action can be 
organized by their UV dimensions. Symanzik improvement then consists of adding 
higher dimension improvement terms to the action, mimicking the effects of the UV 
modes left out on the lattice. To fix the coefficients of these terms, one then proceeds 
as follows. Write down all terms with the appropriate symmetries up to the desired 
order, with arbitrary coefficients. Tune the coefficients by matching to a sufficiently 
large set of observables, calculated either in perturbation theory or non-perturbatively 
in a Monte Carlo simulation. Once the tuning is completed, the coefficients in the action 
will be functions of the physical couplings and a set of redundant couplings (see below) . 
Improving (composite) field operators involves a similar process, which must be repeated 
for each independent field. 

The above program is quite difficult to carry through in practice, not only 
non-perturbatively, but even in perturbation theory. In the past, standard lattice 
perturbation theory suffered from the rather debilitating problem that it did not seem 
to work very well, at least compared to continuum perturbation theory. This has 
now largely been understood [8] as being due to large renormalizations from tadpole 
diagrams, which occur in (naive) lattice but not continuum perturbation theory. Using 
a simple mean-field type method, known as tadpole improvement, one can design more 
continuum-like operators in which the tadpole contribution is greatly reduced. Tadpole 
improvement has been shown to work well for a variety of actions on surprisingly coarse 
lattices [2, 9, 10, 11, 12, 13, 14, 15]. 
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We emphasize that tree-level tadpole improvement should be thought of as a first 
step in a systematic procedure of improving lattice actions. The next step can be further 
perturbative improvement, or, if there are reasons to believe that this is not sufficient, 
non-perturbative improvement. 

It turns out to be substantially harder to improve lattice fermions than gluons, even 
on the classical level. The reason, ultimately, is the first-order nature of the fermion field 
equations, which leads to the well-known doubler problem, which we will discuss later. 
For Dirac fermions (quarks), Wilson [16] solved this problem by adding a second-order 
derivative term to the action. This term breaks chiral symmetry at 0{a). Such errors 
— which are much larger than the 0(0^) errors of "naive lattice fermions" — are too 
large for this action to be useful in coarse lattice simulations. 

The point of this work is to present doubler-free quark actions, for light and heavy 
quarks, that are classically improved to high order. The general tool to construct such 
actions will be field redefinitions; they allow one to introduce second-order derivative 
terms without destroying improvement. To allow the simulation of heavy quarks — 
and also to avoid potential problems due to unphysical branches of the quark dispersion 
relation, which are generic to improved actions — we can use anisotropic lattices [17, 18]. 
Let us discuss these ideas in turn. 

Field redefinitions are just changes of variable in the path integral, so they do 
not affect spectral quantities (at least if one takes into account the Jacobian). Off- 
shell quantities of course do change. Since field redefinitions involve one or more free 
parameters, they lead to so-called redundant couplings, whose values can be adjusted 
at will. This freedom can be used to solve the doubler problem, for example. In 
other situations, in particular on the quantum level, it is very convenient to simplify 
an improved action by setting certain couplings to zero. This leads to the concept of 
on-shell improvement, where only spectral quantities can be obtained directly from the 
action (by improving composite operators one can, however, also obtain their matrix 
elements between physical states). 

The "canonical" procedure for obtaining a doubler-free quark action correct up to, 
say, 0(a") classical errors involves the following three steps:^ 

1. Start with the continuum Dirac action and apply a field redefinition introducing 
even-order derivative terms into the action. 

2. Expand the continuum operators in the transformed action in terms of lattice 
operators up to 0(a") errors; this step will be referred to as the truncation. The 
even-order lattice derivative terms will eliminate the doublers that would be present 
without the field redefinition. One can stop here if one is only interested in spectral 
quantities; they will be classically correct up to ©(a") errors. 

3. To classically also improve off-shell quantities, undo the field transformation (now 
on the level of the lattice action). The resulting action and fields differ only at 

^The simple recipe to follow is a streamlining of what can be found in [6] together with [19]. 
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0(a") form their continuum counterparts, and, in contrast to a naive discretization, 
have no doublers. 

We emphasize that the improved actions so constructed are (classically) improved in 
every respect; the improvement of interactions does not have to be checked separately. 
When applied to lowest order (n = 2) the above procedure gives the Sheikholeslami- 
Wohlert action, originally suggested [6] as an improvement of the Wilson action. The 
next order (n = 4) yields the class of "D234" actions (in addition to the second order 
derivative Wilson term, they also contain third and fourth order derivative terms). 

As alluded to earlier, a problem generic to actions improved beyond 0{a) is the 
existence of unphysical branches of the free dispersion relation, simply due to higher 
order time derivatives in the action. We will refer to these extra branches as lattice 
ghosts. Their energies are at the scale of the (temporal) cutoff, so they will eventually 
decouple as the lattice spacing is decreased. For the lattice spacings used in practice their 
effect on, say, the hadron spectrum has not been thoroughly studied, but they certainly 
affect the renormalization of the improvement terms in the action. In addition, they can 
complicate numerical simulations by introducing oscillations in correlation functions at 
small times. 

If either of these issues turns out to be a problem, one can deal with the ghosts in 
one of two ways. Firstly, one can use field transformations to replace the temporal 
with spatial derivatives. This produces somewhat more complicated actions, as we will 
see. Alternatively and secondly, one therefore might want to use anisotropic lattices 
with smaller temporal than spatial lattice spacing, at < Ug, to push up the energy of the 
ghosts and decouple them. 

Besides effectively solving the potential problem of ghost branches, the use of 
anisotropic lattices has other advantages: 

• By choosing at sufficiently small, one can simulate heavy quarks within a relativistic 
framework [11] without the prohibitive cost of a fine spatial lattice. This provides 
a simple alternative to the NRQCD [20] and Fermilab [21] formalisms. 

• The signal to noise ratio of a correlation function calculated in a Monte Carlo 
simulation decays, generically, exponentially in time. Choosing a smaller at gives 
more time slices with an accurate signal, allowing for more precise and confident 
mass determinations. This is important for particles with bad signal/noise 
properties, like P-state mesons [22] and glueballs [23]. 

• It facilitates thermodynamic studies — one of the reasons being simply that it is 
easier to take independent derivatives with respect to volume and temperature if 
one can vary at independent of Og — especially at high temperatures. 

• It allows for significant simplifications in the design of improved actions. This will 
be relevant for our D234 actions. 
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All these advantages come at a price. Because they have lost part of their axis- 
permutation symmetry, anisotropic actions have more independent coefficients. This is 
not a problem at the classical level, but at the quantum level some of these coefficients 
will have to be tuned to restore space-time exchange symmetry. Large renormalizations 
violating space-time exchange symmetry were in fact seen in first attempts of using 
anisotropic lattices, see [17] and references therein. We find that with an improved 
gluon action and a tadpole improvement prescription appropriate to anisotropic lattices, 
such effects are relatively small, at the level of several percent on coarse lattices [18]. 

So far we have concentrated on classical Symanzik improvement; however, the 
improvement of fermion actions is also more difficult on the quantum level. For pure glue, 
the largest error at order is the violation of rotational invariance, which can be tuned 
to zero non-perturbatively, by demanding rotational invariance of the static potential at 
physical distances.^ Actually, it seems that most of these errors are removed by tadpole 
improvement [2, 18]. 

Wilson- type quarks, on the other hand, have 0{a) errors on the quantum level, and to 
eliminate them one has to tune a term that violates chiral but not rotational symmetry. 
The leading errors behave in the opposite way; they violate rotational symmetry 
but not chiral symmetry (so they are similar to the errors of gluons). The 0{a) and 
(leading) errors of Wilson-type quarks therefore have very distinct effects on spectral 
quantities, and can be tuned iteratively, even on a non-perturbative level, by demanding 
chiral symmetry for the former, and rotational symmetry for the latter. 

As for glue it seems that tadpole improvement does quite a good job in estimating 
the coefficient of the 0{c?) terms that lead to a restoration of rotational symmetry. 
Concerning the restoration of chiral symmetry to eliminate 0{a) quantum errors, Liischer 
et al have recently shown in some beautiful work [7] how to implement this in practice 
for the case of SW quarks on Wilson glue. 

The outline of the remainder of this paper is as follows. In sect. 2 we discuss 
naive lattice fermions, doublers and ghosts. We proceed in sect. 3 to describe in more 
detail the three-step procedure to eliminate doublers, which we then apply to derive 
the Sheikholeslami-Wohlert and D234 actions on a general anisotropic lattice. Several 
special cases and variations are also discussed, including a completely ghost-free D234- 
like action. In sect. 4 we investigate the large mass behavior of the D234 actions. We 
conclude in sect. 5 with a brief summary and sketch of future steps in the improvement 
program. 

Appendices A and B summarize our notation for euclidean continuum and lattice 
QCD, respectively. The reader might want to skim these appendices before starting 

with the main text, and later refer back to them as necessary. Appendix C discusses the 
dispersion relation of the D234 actions. Finally, in appendix D we give some formulas 

^To on-shell improve gluons at order on an isotropic lattice one has to add two terms to the leading 
plaquette, which one can choose to be the "rectangle" and the "parallelogram" [4] . It is a certain linear 
combination of the coefficients of these two terms that can be tuned non-perturbatively by demanding 
rotational invariance. Since the coefficient of the parallelogram seems to be very small (it certainly is to 
one loop) [5, 2], this amounts to an "almost full" non-perturbative tuning of the glue at order . 
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useful in the tadpole improvement of the D234 actions. 

Brief accounts of various parts of this work have appeared earlier in [10, 11, 24, 25]. 



2 Naive Lattice Fermions, Doublers, and Ghosts 

Discussions of Dirac fermions on the lattice^ usually start with the so-called naive lattice 
fermions, specified by the fcrmion operator y + m. Here is the usual first order, anti- 
hermitean, covariant lattice derivative. 



in terms of the link field U^{x) (cf. appendix B for details). This derivative differs at 
order from the continuum Dirac operator = — iA^. 

One way of stating the origin of the doubler problem is that V^^ decouples even and 
odd sites of the lattice. This leads to a doubling (per direction) of the number of degrees 
of freedom on the lattice. If it were not for this problem, naive fermions would provide 
a lattice discretization of Dirac fermions with order errors. Similarly, one could use 
an improved operator, such as 



where is the standard second order lattice derivative of appendix B, and the subscript 
"c" stands for "continuum-like". The fermion operator +m defines what we will refer 
to as naive improved lattice fermions. They would provide a lattice action with only 
order errors — again, if we could ignore the doublers. 

The simplest way of discussing the doubler problem for a generic lattice action is 
in momentum space. Let us consider the dispersion relation E = E[p), obtained from 
the poles of the free euclidean propagator, po = po{p), via E = ±ipo- The two signs 
correspond to particle and anti-particle. For simplicity we will factor out the sign and 
consider only solutions where the (real part of the) energy is positive. The quantitative 
details of the dispersion relations of the actions considered in this paper are discussed 
in appendix C. In figure 1 we show the massless dispersion relations of naive and naive 
improved fermions on an isotropic lattice. There are several noteworthy features. 

• For naive fermions the one branch of the dispersion relation presented in figure 1 
is purely real. Since we can only exhibit a cross section of the energy surface, one 
sees only one of the spatial doublers, which account for half of the doublers. The 

term "spatial doubler" refers to the fact that for each possible energy E there are 
generically eight momenta p (with all components positive) such that E = E{p). 

^It is even more difficult to put chiral fermions on tlie lattice [26]. This we will not attempt. See [27] 
for recent work on this problem. 



Vnij{x) = - — Ui,{x)'tp{x + 1^) - U-i,{x)i;{x - 



(2.1) 



= V^(l - UlA^) 



(2.2) 
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Figure 1: (Real part of the) massless dispersion relation aE = aE{p) as a function of 
a|p|, with p oc (1, 1, 0), for naive fermions (dashed) and naive improved fermions (solid) 
on an isotropic lattice. For comparison we also show the dispersion relation of continuum 
fermions (thin solid). 

• We also can not show that for each possible E = E{p) of naive fermions there 
is another pole of the propagator at E + iir/aQ. This (as well as the existence 
of the spatial doublers) follows from the fact that in momentum space the action 
of naive fermions only depends on af^p^ = sin(a^p^), which is invariant under 
o-i^Pn — ^ TT — a^Pi^. These complex poles constitute the temporal doublers. 

• For naive improved fermions the picture is more complex. There are now four 
branches. The lowest branch is somewhat pathological in that its imaginary part 
is vr/a for all momenta and in that its real part is lower than that of the physical 
branch. It is easy to see that this branch is related to the temporal doubler of the 
naive fermion action. 

Clearly, neither of these two actions corresponds to what one might expect a lattice 
Dirac fermion to look like. Due to the doublers, the naive fermion action actually 
describes 16 Dirac fermion species in the continuum limit [28]. In addition to spatial 
doublers, naive improved fermions have ghosts (or lattice ghosts, if confusion could arise), 
as extra branches of the dispersion relation will be called from now on.^ As mentioned 

*To be sure, these ghost branches are not related to the ghosts appearing in loop diagrams of 
perturbation theory in non-abelian gauge theories. These branches do not describe independent particles; 
they are just lattice artifacts related to the lattice "particle" described by the physical branch of the 
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in the introduction, if the ghosts should turn out to be a problem, there are always ways 
of eliminating them, to be discussed in the next section. 

A solution to the doubler problem was proposed by Wilson [16], who suggested to 
add a second-order derivative term (now known as the Wilson term) to give the fermion 
operator 



on an isotropic lattice, r is the so-called Wilson parameter. It is easily shown (cf. sect. 3) 
that whereas there is one ghost branch for generic r, there is none for r = 1. There is no 
doubler problem for any r > 0. 

The Wilson action with r = 1 therefore solves the doubler problem without 
introducing any ghosts. However, the addition of the Wilson term introduces 0{a) 
errors, which arc too large for this action to be useful on coarse lattices. Sheikholcslami 
and Wohlert [6] described a modification of the Wilson action that has only O(a^) errors 
for on-shell quantities. Their action differs from the Wilson action by a a-F = J2iiu ^t^uFf_iu 
term (cf. appendices A and B for notation), commonly known as the clover term. The 
free dispersion relation of this action is therefore identical to that of the Wilson action. 

In [6] only on-shell improvement was considered. Later it was realized [19] that by 
performing a suitable change of variables on the fields in the Sheikholeslami- Wohlert 
action, one can also calculate off-shell quantities up to C(a^) errors, at tree level. We 
will present a succinct derivation of all this, and its generalization to higher orders of 
improvement on anisotropic lattices, in the next section. 

3 Improved Lattice Fermion Actions 
3.1 Improvement without Doublers 

As seen in the previous section, the naive and naive improved fermion actions have a 
doubler problem. More generally, this is true for any fermion matrix of the form 



which preserves chiral symmetry. (And even more generally, this follows from the 
Nielscn-Ninomiya theorem [26]; see also [28].) To avoid doublers we therefore should, 
following Wilson, introduce chiral symmetry breaking, even-derivative terms (or 
powers thereof) into the action. However, we would like to avoid the 0{a) errors that a 
naive addition of a Wilson term entails. This can be achieved by a field redefinition. 
The simplest way to proceed is to perform the field redefinition in the continuum and 

dispersion relation. What justifies naming them ghosts is that they usually (but not always) give 
negative contributions to the spectral representation of correlation functions. In practice this leads 
to a characteristic "dip" in effective mass plots. 




(2.3) 



+ ...), 



(3.1) 
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only subsequently discretize the action. Starting with the continuum action 

J ^cM^^c = J d^x iPcix) {Ip + mc) Vc(a^) , (3.2) 
we perform a field redefinition 

iPc Mc = i> Mn ^ , Mn = Ct^ . (3.3) 

Note that a field transformation does not affect spectral quantities, at least if we take into 
account the Jacobian of the transformation. Classically the Jacobian does not matter. 
On the quantum level its leading effect is to renormalize the gauge coupling. 
Our canonical choice of field redefinition is (with Qc acting to the right) 

Qc = ^c, ^c^c = . (3.4) 

At this point oq is just a constant with the dimension of length, but in the subsequent 
lattice discretization oq will become the temporal lattice spacing. 
The transformed fermion operator Mq reads 

Mq = rUc + P - ^rao {p^ - ml) 

= mc(l + ^raomc) + p - ^rao (j^ + ^a-F^ , (3.5) 

where we used eq. (A. 11). We can now put the above action on the lattice by discretizing 
]p and -D^ and F^^y to some order, a", using, for example, eqs. (B.9), (B.IO), (B.15) 
and (B.16) in the n = 4 case. Let us call the lattice action so obtained M. 

If one is only interested in spectral quantities, one can use the propagator G = 
in further calculations. Off-shell quantities will then generically have 0{a) errors, since 

= 1 + 0{a) and therefore = ipc + 0{a). However, as our third step, we can also 
improve off-shell quantities by undoing the field transformation. To do so, we use the 
obvious lattice versions of the operators in eq. (3.4), which differ from them at order 
a". Let us call these operators Q and Q. The action obtained by undoing the change of 
variable is M using fields that differ from the original continuum fields only at 
0{aP). The propagator of this action is 

G = n M-^ n = M-^ + 0(o") . (3.6) 

Note that undoing the field transformation on the lattice does not lead to the 
(re) appearance of doublers. 
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Figure 2: The energy asE{p) of free SW/Wilson fermions with r = 1 as a function of a^lpl 
with p oc (1,1,0). On the left wc show the massless case on 1:1 (sohd), 2:1 (dashed), 
and 5:1 (dot-dashed) lattices, as well as continuum fermions (thin solid). On the right 
we show the same for mass agnic = 1. 

3.2 The Sheikholeslami-Wohlert Action and 0{a) Terms 

Using the leading discretization of the derivatives in (3.5) gives 

Msw = mc{l + Iraorric) + f- ^rao - ^raoa-F . (3.7) 

This is the Sheikholeslami-Wohlert (SW) action on an anisotropic lattice. For F^j, one 
can use the so-called c/ower representation (cf. appendix B), which has O(a^) errors. By 
construction this action has classical O(a^) errors for spectral quantities; also for off-shell 
quantities if we undo the change of variable. To obtain the Wilson action one must by 
hand set the clover term, a-F, to zero, thereby incuring an O{ao) error in the presence 
of a non-trivial gauge field. 

This action has no doublers for any r > 0; it has no ghost branches either if and only 
if r = l, which is therefore the canonical choice. For future reference we show in figure 2 
the free dispersion relations for the SW/Wilson action with r = l for various anisotropies 
Us /at and masses. 

Recall that in general more operators are needed for quantum improvement (even 
on-shell) than for classical improvement. For Wilson- type quark actions on an isotropic 
lattice, however, it is easy to see that the clover and Wilson terms are the only ones 
allowed at 0{a) by gauge and (discrete) rotational symmetry [6]. Their coefficients of 
course renormalize on the quantum level. By a field transformation of the canonical form 
discussed above one can adjust the coefficient of one of these terms to any desired value. 
It is natural do so for the Wilson term; in the SW case to maintain its "bare" coefficient 
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at the canonical value r = 1, for example. To eliminate quantum 0{a) errors one then 
has to tune the coefficient of the clover term. 

On an anisotropic lattice the situation is slightly more complicated. The allowed 
operators at 0(a) consist of the spatial and temporal parts of the Wilson and 
clover terms, and the additional operator [yoDo,J2i"fiDi]. The most general field 
transformations fi,, and allowed in this situation lead to three redundant operators, 
so that one has to tunc two coefficients at 0{a). These can be chooscn to be the spatial 
and temporal parts of the clover term. Note that on an anisotropic lattice one must 
also allow a relative coefficient between the temporal and spatial kinetic terms at 0(a°), 
which can be tuned non-perturbatively by demanding a relativistic dispersion relation 
for the pion, say, at small masses and momenta. 



3.3 The D234 Actions 

Going to the next order in the expansion of the continuum derivatives in (3.5) gives the 
class of D234 actions 

Md234 = rnc{l + ^raomc) + ^ T^^V^ (1 - fe^^a^A^) 

- I rao + i ct-f) + ^ c^a^ (3.8) 

where, at this point, 

^'^ = 6' = ^ ^'-'^ 

The specific D234 action defined by the coefficients in (3.9) will be referred to as 
"D234c(r)", where "c" refers to the fact that this action is obtained by our "canonical" 
field redefinition without any further modifications. If we use an improved representation 
of the field strength, as in eq. (B.16), this action only has 0{a^) classical errors. There 
is no canonical choice of r for this action. It will generically have three ghost branches, 
as illustrated in fig. 3 for the case of r = 1 (for r = 2 there arc only two ghosts, but the 
lowest one is too low for this choice to be interesting, except perhaps on very anisotropic 
lattices). 

It remains to be investigated what effect these ghost branches have on the quantum 
level, but one should certainly consider designing actions with fewer and higher-lying 
ghost branches. As we will see, the necessary "tuning" of the D234 actions introduces 
classical errors in addition to the 0{a^) ones. However, this is probably irrelevant, since 
the latter errors are unlikely to dominate on the quantum level anyhow. 



3.4 Tuning the D234 Actions 

To investigate the ghosts let us study the free dispersion relation corresponding to the 
general D234 action (3.8). The details are discussed in appendix C, to which we refer 
for the proof of any non-obvious facts we will use. For generic r, and the dispersion 
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Figure 3: As in figure 2, for the massless D234c(l) action on 1:1 (sofid) and 2:1 (dashed) 
lattices. We only show the real part of the energy, relevant for the top branch of each 
anisotropy, which has imaginary part tt/oo, and to the right of the branch point on the 
1:1 lattice. 

relation is a quartic equation for sinh(ao£'/2)^, E = E(p), so there will be three ghost 
branches.^ Since the coefficients of the quartic are real, the energies will be real, come in 
complex conjugate pairs, or have imaginary part ir/aQ. Note that the qualitative branch 
structure {e.g. the number of branches) depends only on the temporal coefficients r,bo 
and cq. For m,c = and p = the only way a lattice spacing enters the dispersion relation 
is via uqE. This implies that for small momenta and masses the height of the ghosts is 
inversely proportional to the temporal lattice spacing. 

The most basic question we can ask about the ghosts, is how many ghost branches 
we can completely eliminate by a suitable choice of the free parameters. We summarize 
the conclusions of appendix C concerning this question as follows: 

• If we choose bo = 2co there will be at most two ghosts. 

• If we further choose r = 1 — 2bo or 6o = there is (at most) one ghost. 

• The only way to eliminate all ghosts is to choose r = l,bQ = cq = 0, which is of 
course the standard SW/ Wilson case (if r 7^ 1 the SW/ Wilson action will have one 
ghost). 

^Remember that we do not count the particle anti-particle symmetry E <-^ —E in the number of 
solutions. 
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If we want to go beyond the SW action we will therefore have at least one ghost 
branch. Let us now discuss our "favorite" D234 actions on isotropic and anisotropic 
lattices in turn. 

3.4.1 Isotropic D234 Action 

On an isotropic lattice one will presumably prefer an isotropic action with manifest 
space-time exchange symmetry, so that one does not have to restore this symmetry by 
tuning on the quantum level. With = ^, the above results imply that when we require 
one ghost we must choose r = | and = We will refer to this action as the "isotropic 
D234" action. It was introduced in [10]. 




Figure 4: As in figure 3, for the massless isotropic D234 (solid) and SW (dotted) actions. 

Since the violate eq. (3.9) this action has classical 0{a^) errors. At zero mass and 
momentum the one ghost branch of this action is at aE{Qi) =ln 7 1.9459. Its massless 
dispersion relation is shown in figure 4, together with that of the SW/Wilson action. 

We will see in sect. 3.5.3 that, like the SW and Wilson actions, this action can be 
coded very efficiently using the "projection trick" . 

3.4.2 Anisotropic D234 Action 

On an anisotropic lattice we can have one ghost with only 0(aQ,a^) errors simply by 
modifying the coefficient cq of the D234c(|) action to he cq = We will refer to this 
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action by the name D234i(|).^ Its dispersion relation on a 2:1 lattice is compared in 
figure 5 with that of the SW action. Note in particular the impressive dispersion relation 
in the massive case, indicating that such an action might be very useful for heavy quark 
simulations. 




Figure 5: As in fig. 2, for D234i(|) on a 2:1 (solid) and SW/Wilson fermions on a 1:1 
(dotted) lattice. The real part of E{p) is shown to the right of the D234i(|) branch 
point. In the massive case the bare masses are tuned so that asE{0) = l. 

We should point out that the restriction of the D234i(|) action to an isotropic lattice 
does not give the isotropic D234 action. Because the spatial Cj of the former were chosen 
to not introduce af errors, this action has anisotropic coefficients even on an isotropic 
lattice; it was designed for use on anisotropic lattices. 

3.4.3 Variations 

By relaxing the requirement of just one ghost one can construct actions that might be 

interesting for sufficiently anisotropic lattices. We will not discuss these here, but just 
make the general remark that for larger anisotropics as /at it is advantageous to choose 
larger values of r. Otherwise one will recover spatial doublers, as is obvious from the 
fact that our canonical field transformation (3.4) J7c ^ 1 as oq ^ for fixed rJ 

We should briefly discuss one modification of the D234i(|) action. Namely, it 
is possible to construct a very similar D234 action that in the free case has only 
0(0^) classical errors, with ©(oq) errors entering only in the presence of a gauge field. 

®The "i" stands for mass- "independent" , since the coefficients in this action enjoy this property. In 
the next subsection we will describe a closely related D234 action with mass-dependent coefficients, to 
which we have previously [11] given the name D234(|). 

'^If one avoids the reappearance of doublers by letting rao approach some non-zero limit as oo — » 0, 
one will obtain one ghost branch, with energy 2/rao for small masses and momenta. 
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This action is obtained by a somewhat more complicated change of variable from the 
continuum Dirac action, 

dc = n,, = 1- Kao {p-mc- K,) , K, = J2 Cm^mTm^m (3-10) 

with free parameters in addition to r. Requiring one ghost and O^Gq) errors only in 
the presence of a gauge field leads uniquely to the choice of coefficients 

1 1 

6 1 + faomc 

1 1 + laonic 
6 1 + ^aorric 

I. 

2 1 + faomo _ ,3 

3 1 + j^aorric 

as well as = 0, 6j = ^ and q = rao/24aj. Wc will refer to this action as D234(|) 
(labelling it by the value of r at mc = 0). This was the action used in the anisotropic 
lattice simulations described in [11]. Note that numerically this action differs significantly 
from D234i(|) only for very large masses; for mc = they are identical. 

Finally, we remark that it is easy to improve the dispersion relation of D234-like 
actions at large momenta still further by introducing suitable fifth and sixth order 
derivative terms. This can be done by a field transformation and/or by simply adding 
such terms to the action. However, in the relevant momentum regime the hadron 
dispersion relations measured in simulations of D234 actions are already so good, even 
on coarse lattices [10, 11], that the additional cost and complication from the higher 
derivative terms seems unjustified. 

3.5 Other Actions 

We conclude this section by briefly discussing several other classes of actions (or other 
ways of writing actions) that are of interest for various conceptual and practical reasons. 

3.5.1 Ghost-free D234-like Actions 

We will now demonstrate that it is straightforward to write down a highly improved 
action, at tree level, that has no ghosts whatsoever. This comes at a price, of course. 
Such an action is more complicated and therefore more costly to simulate. 

The idea is to use field transformations to eliminate the cubic temporal terms 70V0 Aq 
in the naive improved fermion action in favor of spatial terms. Starting with the 
continuum action Mc = rric+p, we perform a field transformation ipc = ^ci V'l V'c = V' ^ci 



eo = 

bo = 
Co = 
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with^ 



Oci = 1 + 



n. 



cl 



1 + 



12 
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Dl - {p - mc)iV> + rric) 
D^^-{y) + mc)i]p-mc) 



where the purely spatial derivative y> = J2i JiDi. This gives 



(3.12) 



(3.13) 



Note that, when discretized, the term p + ^Oq^oDq will not contain any lattice time 
derivative above the first, and there is no other temporal derivative in the action. We 
can now proceed with a second change of variable to introduce even derivatives, defined 
by ^c2 = ^c2 and 



^c2 = 1 + |rao + ^ag7o£>o - "^c - SKc) . (3.14) 



This implies 



fic2 ^ci Mc r^ci ^c2 = rnc{l + hraomc) + Jp + -aoTo-Dg + (1 + raouic) 6K, 
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Finally, we discretize this action. Using 



(3.15) 



(3.16) 



we find 



MD234gf = mdl + Iraonic) + + Yc + (1 + raonic) 5K 

- \ra, (E + k • ^) + 1^ E ' (3-17) 

M i 

SK = ^{ml-Y^Ai-^a-F,Y + mc} = 6K, + 0{a^), (3.18) 



where = 70 Vo is the unimproved temporal lattice derivative, is the improved 
spatial derivative (see (B.12)) and for the field strength F one should use an improved 
discretization. This action has O(ao,a^) errors in on-shell quantities. As usual, by 

*This field transformation is similar to ones used in [21]. 
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Figure 6: As in figure 2, for the ghost-free D234-hke action on 1:1 (solid), 2:1 (dashed) 
and 3:1 (dot-dashed) lattices. 

undoing the above field transformations one can achieve the same errors for off-shell 
quantities. 

Obviously the above action is the analog of the D234 actions of the previous sections, 
with SK appearing in place of cubic and quartic temporal derivative terms. Now, 
however, we have just one branch if we set r = 1, as for the Wilson and SW actions. 
This is illustrated in figure 6. 

For the aj] errors to be negligible compared to the ones, we can again use anisotropic 
lattices. In that case the ghost branches of the D234 actions of the previous section 
are presumably harmless, and it seems doubtful that having no ghosts outweighs the 
disadvantage of having to include the costly anti-commutator term 6K. Comparisons 
to simulations with this action should, however, allow one to discern whether the ghost 
branches have any effect besides that on correlation functions at small times. 

Another, less ambitious ghost-free action can easily be constructed if one is willing to 
tolerate Og in addition to Og and errors. Such an action can be obtained, for example, 
by simply neglecting the temporal third and fourth order terms in the D234c(l) action. 

3.5.2 D34 Action 

Recently a tadpole-improved version of an improved action discussed in [29] was used in 
a Monte Carlo simulation of the hadron spectrum [30]. This action has third and fourth 
order derivative terms, but no second order Wilson term. We will therefore refer to it as 
the D34 action. Generalized to an anisotropic lattice it reads in our notation 

Md34 = mc + % ■ (3-19) 
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Figure 7: As in figure 5, for D34 fermions with = | on 1:1 (solid) and 2:1 (dashed) 
lattices. Note that this action has four branches. They consist of two pairs of complex 
conjugate solutions for sufficiently large momenta. For small momenta there are usually 
two real branches. For asmc = l on an isotropic lattice, however, all branches axe complex 
even at small momenta. 

In [30] an isotropic lattice with all = g was used. The D34 action is obtained from 
the naive improved fermion action not by a change of variable, but simply by adding the 
term, which leads to an error (and means that the fields are already improved up 
to errors; no change of variable has to be undone). This action is of course a special 
case of the general class of D234 actions, with r = and 6^ = g. It has three ghosts, 
except for cq = when there are two. There is no obvious canonical choice for the c^; 
for all values of the ghost branches seem to lie rather low. In figure 7 we show various 
dispersion relations for the = g case. 

3.5.3 W-Actions 

Recall that the Wilson and SW actions with r = 1 can be efficiently coded using the 
"projection trick" , which exploits the fact that the "Wilson operator" [31] 

= Vf,- ^7mAm (3.20) 

can be expressed in terms of (spinor-) projection operators. Most conveniently this is 
expressed as 

li.W^ = -V+P- + V;P+, (3.21) 

where Pj^ = ^(1±7^) are projectors on two-dimensional, orthogonal subspaces (for fixed 
/J,), and 

V±^(x) = ±— (c/±^(x)V(x±/x)-V(x)) (3.22) 
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are forward and backward derivatives. 

Actions that can be expressed as (low-order) polynomials in the W^, 

M = Y.— Q^ia^i^W^) (3.23) 

(with the possible addition of mass and clover terms, cf. below) will be cheap to code. 
Another advantage of such an action is that one can immediately read off the number 
of branches of the dispersion relation: If Qq{x) = x + . . . is an n-th order polynomial, 
there are exactly n branches, counting possible degeneracies (the proof of this statement 
is left as an exercise to the reader). Note that if all Q^(x) are n-th order polynomials, 
the class of such "W23. . . ra" actions is a proper subset of the class of "D23. . . 2n" actions 
considered previously, which generically have 2n branches. 

It is interesting to ask when such W-actions are improved. We first remark that the 
continuum derivative can be written as 

a^,l^,D^, = - ln(l - 0^7^!^^) = a^-i^M^^ + \{a^,-f^W^,f + ^{a^,-i^W^f + . . . (3.24) 

Truncation of this expansion does, however, not seem to lead to promising (on- and off- 
shell) improved actions, since the expansion converges too slowly and the ghost branches 
lie too low. 

Alternatively, we may ask when a W-action is only on-shell improved, related to 
some order in a by our canonical field transformation (3.4) to the continuum Dirac 
action. Restricting ourselves to isotropic actions from now on - the improved actions 
in question are not simply expressed in terms of VF^ on anisotropic lattices — we write 

M = mc(l + \ram^) + - ^ QiaiuW^:) - \ raa-F . (3.25) 

For a first-order polynomial, Q{x) = x, wc recover, of course, the SW action for r = l, 
with O(a^) errors. To see which higher order polynomials correspond to improved 
actions, we expand out the in terms of and and compare with the results 
of the previous subsections (specifically, the D234c(r) actions and their higher order 
analogs). 

The second order polynomial Q{x) = x + ^x^ with r = ^ gives a "W2" action, 
which has O(a^) errors. It is identical to the isotropic D234 action [10] of sect. 3.4.1. 
This implies that the application of the isotropic D234 operator is only about twice as 
expensive as that of the r = 1 SW or Wilson operator. 

For a third order polynomial one can reduce the errors to O(a^) by choosing 
Q{x) = X + ^x^ + ^x^ and r = to give a "W23" action (one must now also use an 
improved F in the clover term). This action is equal to a D23456 action with r = ^, 
ft/i = g) C/j- = dfj- = ^5 e^t = ^) where d/j, and are the coefficients of fifth, 
respectively, sixth order derivative terms, defined analogously to and c^. This action 
has two ghost branches; for all masses and momenta they form a complex-conjugate pair. 
For zero mass and momentum their energies are aE{0) 1.528 ± 0.897 z. 
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4 Large Mass Behavior of the D234 Actions 



Our D234 actions were initially designed with light quarks in mind. In contrast to [21], 
for example, where renormalization conditions are imposed on mass-shell, even for heavy 
quarks, we effectively do so at zero quark mass, since we always expand in powers of 
a^, assuming a^p^ to be small. Heavy onium systems are non-relativistic, so masses are 
large but momenta arc small. One would therefore expect that by going to an anisotropic 
lattice, where atnic is small, our D234 actions could be used to simulate such systems 
on lattices where the spatial lattice spacing is still quite large. Decreasing the temporal 
lattice spacing by a factor of, say, 3 — 5, would present a relatively minor increase 
in cost, given the exciting prospect of obtaining accurate results for e.g. charmonium 
within a relativistic framework. The charmonium system is difficult to simulate. It is 
light enough for the NRQCD expansion to become problematic, but too heavy for the 
usual light quark actions on isotropic lattices to be accurate. A highly improved quark 
action on an anisotropic lattice seems taylor-madc to simulate such a system. 

To get a quantitative idea of how small we have to choose at, let us investigate in the 
free case when the D234 actions break down at large masses. As indicators of break-down 
we will consider E{0)/mc and the "effective velocity of light" c(p), defined by 

P 

where E{p) is meant to be the physical branch of the energy. We know that for small 
masses E{0)/mc and c(0) are 1 up to order O(afm^) or O^afmf) corrections. For large 
masses, it is clear from the figures presented in sect. 3, that for the relevant D234 
actions the two lowest real branches of the energy (at fixed momentum) will eventually 
merge. (For zero momentum this branch point occurs at the same value of atruc for any 
anisotropy. For non-zero momentum this becomes true only asymptotically, for large 
anisotropy.) At the branch point c(p) diverges. 

In figure 8 we compare c(p) for the D234i(|) and SW/Wilson fermions for various 
anisotropies. We show c(p) as a function of mass for both p = and the (quite large) 
value p = (1,0, 0)/as. 

As expected, c(p) stays close to 1 for a larger and larger mass range as the anisotropy 
increases. Note how much better the D234 action behaves in this respect than the 
SW/Wilson action, in particular at non-zero momentum. It is also interesting to observe 
that for the SW/Wilson action c(p) decreases more or less monotonically as a function 
of both nic and |p|, whereas for the D234 action c(p) increases as a function of nic for 
fixed p, but (more or less) decreases as a function of |p| for fixed nic (an isotropic lattice 
provides an exception to this latter statement). This is presumably one of the reasons 
we find in Monte Carlo simulations that the D234 actions exhibit an excellent dispersion 
relation for various hadrons up to surprisingly large masses and momenta. 

As is clear from figure 8, for the D234 actions the transition from c(p) ~ 1 to 
divergence is quite rapid. This provides a sensitive indicator of the breakdown of the 
D234 actions at large masses. When c(p) begins to get larger than 1, the D234 actions 
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Figure 8: c(p = 0) (left) and c{p = {l,0,0)/as) (right) as functions of a^mc for D234i(|) 
(solid) and SW/Wilson (dotted) fermions on 1:1, 2:1 and 3:1 lattices. 

should not be trusted anymore. When this happens, it should only be necessary to 
decrease at by a relatively modest amount to be able to simulate the quark mass of 
interest. 

To reinforce this last point, let us mention that among the anisotropic D234 actions 
we discussed, the finite-mass errors are largest for the D234i(|) action, where 



For comparison, the corresponding expansions for the free SW/Wilson action are 



We see that not only do the D234 actions have a larger power in their scaling errors, 
they are also blessed with small coefficients. The coefficients in the above expansions for 
various other actions can be found in appendix C. There we also discuss the error that 
arises in the presence of a non-zero gauge field. Again we find that the coefficients of the 
error terms are very small as long as atruc < 0{l). 

For the large mass error in c(0) to stay below the one or two percent level, one must 
choose atiric < 0.5 for the D234 actions. For the SW action one must satisfy the much 
tighter constraint atrUc < 0.2 to achieve the same error. Given this bound for the D234 




(4.2) 



(4.3) 



E{0) 




c(0) 



(4.4) 
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actions, we can attain a meson mass 

mqq 2mc = — atf^c (4-5) 

in the charmonium range {mqq 3.0 GeV) with a spatial lattice spacing of !=a 
600 — 1000 MeV and anisotropy = 3 — 5. For the SW action the lattice spacings would 
have to be almost three times smaller with the same anisotropy, making a Monte Carlo 
simulation roughly two orders of magnitude more expensive. Needless to say, on an 
isotropic lattice a simulation of charmonium in a relativistic formalism would be more 
expensive by many, many additional orders of magnitude. 

5 Discussion and Conclusions 

We have seen that using field transformations it is easy to design highly improved, 
doubler-free quark actions at tree level. We emphasize that the improvement includes all 
interactions between quarks and gluons, not just the free quark dispersion relation. We 
have discussed several actions, representing different compromises between the conflicting 
aims of high level of improvement, absence of ghost branches, and simplicity. The 
actions we are currently using in simulations are the isotropic D234 (sect. 3.4.1) and 
the anisotropic D234i(|) (sect 3.4.2) actions. However, given the ingredients provided, 
the reader can concoct many more flavors of actions. 

The actions we constructed arc jiist moderately more expensive to simulate than the 
SW action. For example, the isotropic D234 action is only about twice, the generic D234 
action about four times as expensive than the isotropic SW action with r = 1 (a slight 
additional overhead might be incured in preconditioning). 

The next task, of course, is to put some quantum flesh on these classical bones. 
Exploratory quenched simulation results using tadpole improvement have appeared [10, 
11], and are very encouraging. Further work is in progress [22]. In terms of the general 
improvement program, the next step is to check if a non-perturbative tuning of the 
low-order coefficients in the action is necessary. To eliminate all errors up to 0{a), one 
would only have to tune the clover term, for an isotropic action. It was recently shown [7] 
how to implement this in practice, by demanding the restoration of chiral symmetry at 
zero quark mass. In [7] the case of SW quarks on Wilson glue was considered, but the 
methods apply equally well to other actions. 

As mentioned in the introduction (cf. sect. 3.2 for more details), for anisotropic 
actions there arc more coefficients that can renormalize on the quantum level, and 
in principle have to be tuned in order to restore space-time exchange symmetry. We 
hope that suitable tadpole improvement on anisotropic lattices [18] will reduce these 
renormalizations so that their effect on physical quantities is only on the few percent 
level, where they can be neglected; but this is an issue to be decided by empirical study. 

Already at 0{a^) one might have to introduce a "bare velocity of light" into 
the action, to restore space-time exchange symmetry at leading order. The required 
renormalization is easily determined non-perturbatively in this case, by measuring 
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the dispersion relation of a pion, say, at small masses and momenta. For on-shell 

improvement at 0{a) there is one additional coefficient compared to the isotropic case, 
which one can choose to be the relative coefficient of the temporal and spatial parts of 
the clover term. Using the methods of [7] with a suitable background field, for example, 
it should be possible to determine this renormalization. Finally, we should point out 
that on an anisotropic lattice the effect of these renormalizations is suppressed, since the 
overall coefficient of the 0{a) terms is smaller by a factor ot/a compared to the isotropic 
case (for fixed r). 

In summary, the D234 actions show great promise for accurate QCD simulations 

on coarse lattices of improved glue, for both light and heavy quarks. The low-order 
terms in the actions require further study, and, if necessary, should be determined non- 
perturbatively. We hope that, as for improved glue, quark actions with negligible 0{a) 
and only small quantum O(a^) errors will give accurate results for all observables on 
coarse lattices. The next task of our program will be the determination of renormalization 
constants like Za and Zy for the D234 actions. Once quenched QCD is satisfactorily 
understood, it is time for simulations of full QCD. Realistic simulations of full QCD 
should finally be possible on coarse lattices. 
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Appendix A Euclidean Continuum QCD 

We write the action of euclidean SU(A^) gauge theory in four dimensions as 

Sg[A] = -L y d^a; TV F^,(x)F^,(x) , (A.l) 

where F^,^{x) is the su(A^)-valued hermitean field strength. In our conventions the 
covariant derivative is written in terms of the hermitean gauge field An{x) as 

= d^-iA^, (A.2) 

so that 

F^^ = i[D^,D^] = d^A^-d^A^-i[A^,A^] . (A.3) 

In terms of traceless hermitean su(iV) generators Ta, a = 1, . . . , A'^^ — 1, normalized by 
Tr{TaTb) = ^Sab, we write 

Fi^Ax) = F^Ax)Ta, A^{x) = Al{x)Ta . (A.4) 

The parallel transporter from a point x to y along a curve Cyx is 

U[Cyx] = Pexpfi I' dx'^A^{x')\ gSU(A), (A.5) 

where P denotes path ordering, stipulating that in a series or product expansion of the 

exponential the fields at a point earlier on the curve are to be placed to the right of fields 
at later points. Under a local gauge transformation, A(.) G SU(A''), 

) ^ A(x)F^^(x)A-^(x) 
lJ\Cyx\ ^ k{y)U\Cyx\t^-^{x) . (A.6) 

Recall that the inhomogenous transformation law of A^{x) has been designed so that 
the parallel transporter transforms covariantly as above. 

The action of a Dirac fermion coupled to a gauge field is 

Sf[ip,ip] = J d^x i){x) {Ip + m) ipix) , (A.7) 
where p = J^/j, Ifi^iJ. terms of the euclidean gamma matrices defined by 

7/. = tJ 

{7;., 7.} = 2,5^. • (A.8) 

The spinor fields jp{x) and -ip also carry a suppressed color index in the vector 
representation of SU{N). Under a gauge transformation 

ip{x) —I- A{x)ip{x) 

i;ix) i;{x)A~^ix) (A.9) 
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and Dn'tp{x) transforms like ^'(a;) (again, by construction of An{x)). 

We will need some identities involving the hermitean matrices defined by 

cr^tu = - 2 [Tm' 7i/] = -(^i^fi , or 7^7,, = + ia^^ . (A. 10) 

This implies 

^{lnD^,-f^D^} = 6i,^Dl + -a^^Fi,^ (A.ll) 
Appendix B Lattice QCD 

We will consider a four-dimensional euclidean hypercubic lattice of extent and lattice 
spacing in direction /j, = 0,1, 2, 3. Since we will always work in euclidean space 
(except when we consider dispersion relations) the use of /x = to denote the euclidean 
time direction should not cause confusion. 

Points are labelled by x, y, . . ., as in the continuum. When all spatial lattice spacings 
are identical, they are denoted by a^; the temporal lattice spacing will then be denoted 
by at = ao. We will refer to ^ = ag/at as the anisotropy of the lattice, and sometimes call 
such a lattice a : 1 lattice". For an isotropic lattice we set a = a^. 

Wc use for spatial indices, and boldface letters for spatial vectors. The 

notation x ± is a shorthand for x ± a^/t, where /t is a unit vector in the positive 
/x-direction. 

When working on anisotropic lattices one should in principle be very careful in 
specifying the lattice spacing errors of various quantities. To avoid cumbersome notation, 
we will be careful only to distinguish errors of the form ©(a^) from 0(a") errors, where 
the latter denotes any errors that are not of the form 0{a^) or are a sum of such terms 
with different /x. 

Gauge field dynamics on a lattice is expressed in terms of the link field U^{x), which 
takes values in the gauge group. Un{x) is a parallel transporter from x + /j, to x, along a 
straight line. In terms of an underlying continuum gauge field An{x) we therefore have 
(cf. appendix A) 

U^{x) = fj^ dx'^A^ix')^ , (B.l) 

where now (because we changed the sign of the path compared to eq. (A. 5)) fields at 
points later on the path are to be placed to the right of earlier fields. We will employ 
the notation U-^{x) = Ufj^ix — /j,)^ for the parallel transporter from x — /i to x. 

Using the link field U^{x) is the only known way of maintaining exact gauge invariance 
on the lattice (not preserving manifest gauge invariance would necessitate a costly tuning 
of the various gauge couplings). Under a gauge transformation 

U±^{x) ^ K{x)U^,{x)K-\x±ii). (B.2) 
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With the help of the link field it is trival to construct gauge-covariant first- and 
second-order lattice derivatives via 



1 
2^ 



U^(x)'4>{x + n) - U-^(x)i;{x - ii) 
U^{x)4^{x + /x) + U-^{x)^{x - /x) - 2i^{x) 



(B.3) 
(B.4) 



The operators and are the building blocks of our lattice fermion actions. 

It is useful to observe that in momentum space (with = 1) they correspond to 



A,, 



-pI , aixPix = 2 sin(a^p^/2) 



(B.5) 
(B.6) 



Momenta on a lattice can be taken to lie in the Brillouin zone, defined by — 7r/a^ < P/x < 
vr/o^. Note that p^, in contrast to p^, has a "doubler problem", i.e. it vanishes at the 
edge of the Brillouin zone a^jo^ = vr. 

Since V^^ and A^ are gauge-covariant, many identities involving these operators are 
easy to prove by going to a gauge in which 11^ = 1, where such identities reduce to 
relations between p^ and p^. For example, 



and 



Other identities we will need follow from 

„2 _4 _6 

P^^ = P^.[l + -fP^. + tt^P^ + t^P^ + 



a? 



90 



560 



The operator V^^ has O(a^) errors compared to the continuum derivative. 



(B.7) 
(B.8) 

(B.9) 
(B.IO) 

(B.ll) 



Eq. (B.9) implies that a more continuum-like first-order covariant derivative, involving 
only next-nearest neighbor sites, can be defined as follows 



(B.12) 



The clover operator is a well-known lattice representation of the field strength that 

agrees with the continuum F^y up to 0{a?') errors. To define it, let us expand on the 
path notation used for the parallel transporter in appendix A. Namely, let us introduce 
a shorthand for lattice paths, defined recursively so that Cyxin) denotes the path from 
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X + fj, to X to y, and x = Cxx- Path ordering implies that C/[Cyj;(/i)] = U[Cyx]Un{x). Of 
course, U[Cxx] = ^- 

The clover operator we use, fI}^\x), can be expressed in terms of the sum of the link 
fields around the four plaquettes bordering on x (all counter-clockwise, say), as 



1 



T U[x{f,)iu){-^^){-l^)] + C/[xH(-/x)(-i.)(/x)] + 



where T{M) is the (color-)traceless imaginary part of an AT x AT matrix, 

r(M) = ^{M-M'') - -^ImTrM. 
A simple calculation shows that 



F^u{x) + ^{alDl + alDl)F^,{x) + 0{a' 



(B.13) 



(B.14) 



(B.15) 



Using a suitable discretization of the second-order derivatives acting on the field strength, 
we can therefore define a more continuum-like lattice field strength via 

Fe^,(x) = ^-Fj;^){x) -^-\u^{x)Fj;^\x + tx)U.^{x + tx) + 



u^^{x)Fff{x - ix)u^{x -^J)-{^,^y) 



(B.16) 



Appendix C Dispersion Relation of the D234 Actions 

In terms of p and p introduced in appendix B let us define 

m{p) = mo + -raoY,Pl + Y.^l^"'lpt- 
mo = mc(l + ^raouic) ■ (C.l) 

The inverse of the free propagator of the general D234 action (3.8) is then simply 
ip + rh{p) in momentum space, implying the dispersion relation 

p2 + m(p)2 = . (C.2) 

As a dimensionless measure of the energy it is useful to introduce 

y ^ y{E) = -alpl = Asmh\^) , (C.3) 
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so that 



2-2 ri , y\ 

-alpl = y{l + l){l-boyf . (C.4) 
Expressed in terms of y the dispersion relation for the energy E = E{p) therefore reads 



where 



y (1 + |) (1 - bo yf = alp' + (/x(p) - ^ y + cq y'Y , (C.5) 



/x(p) = aomo + ^ral^pf + ao X] Cjofpf • (C.6) 



Note that the use of y factors out the particle anti-particle symmetry E ^ —E, since y 
is invariant under it. The dispersion relation is generically a quartic equation for y with 
real coefficients. This implies that generically there will be four solutions for ±E, which, 
at a given p, are either real, come in complex conjugate pairs, or have imaginary part 
ir/ao. 

Let us now discuss when the quartic reduces to some lower order polynomial. We 
see immediately that the quartic term cancels when cq = 560. Assuming this holds, the 
cubic term cancels if 6o = or r = 1 — 26o- The remaining quadratic equation reads 



-y' 



60(2 + /x(p)) + \{r'- 1)] +y[l + rfi{p)] = nipf + alp' . (C.7) 



Note that the only way to obtain just one branch is to choose r = l,bo=co = 0, which 
corresponds to the Wilson or SW action. 

Modulo the slight changes necessary for the ghost-free D234gf action, all dispersion- 
related results and plots presented throughout this paper can be obtained by analytical 
or numerical manipulation of eq. (C.5). For completeness we include some more 
results concerning the small mass expansions of E{0)/mc and c(p) defined by -E'(p)^ = 
E{0)' -|-p^c(p)^. Here E{p) is the physical branch of the dispersion relation (cf. sect. 4); 
c(p) is defined only for (sufficiently small) momenta and masses, where there is a well- 
defined physical branch. 

The first few coefficients in the expansions 



m 



rrir 



1 + J2 Ena^m^ 



n=l 

00 



c(0) = 1 + E Cna^m^ (C.8) 



n=l 



are presented in tables 1 and 2 for various actions discussed in sect. 3. These results are 
easily derived by an iterative series expansion of eq. (C.5), using a symbolic manipulation 
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Action 


El 


E2 


E3 


Ea 


E5 


Eq 


SW 

D34(co = i) 
D234 = W2 
D234i(|) 
D234(§) 

D234c(|) 










1 

6 








1 

8 

1 

6 

1 

18 
1 

18 





1 

20 

1 

30 

1 
270 

1 
270 
1 

30 
1 

30 




1 

36 

5 
324 

5 
324 

5 
216 
1 

54 


1 

56 

29 
252 

41 
6804 

41 
6804 

115 
6048 

37 
2268 



Table 1: Expansion coefficients of E{0)/mc, eq. (C.8), for various actions. 



Action 


Ci 


C2 


C3 


C4 


C5 




SW 





1 

3 


1 

4 


7 
180 


5 

48 


1009 
7560 


D34(co = |) 








1 

3 


1 

10 


1 

12 


22 
63 


D234 = W2 








1 

9 


7 
270 


7 
324 


131 
3402 


D234i(|) 








1 

9 


7 
270 


7 
324 


131 
3402 


D234(|) 











1 

10 


5 

72 


1135 
18144 


D234c(|) 











1 

10 


1 

18 


10 
189 



Table 2: Expansion coefficients of c(0), eq. (C.8), for various actions. 



program. Note that E{0)/mc and c(0) are independent of the spatial A? terms. For the 
D34 action we therefore only have to specify cq in tables 1 and 2. This also implies 
that D234i(|) and the isotropic D234 action have the same E{0)/mc and c(0) expansion 
coefficients, cf. tables 1 and 2. 

For the SW action with r = l the exact mass E(0) is given by the simple formula 

ao-E(O) = In 1 + aomdl + ^aQUic) ■ (C-9) 

We leave it as an exercise to the reader to derive analytical results for various other 
special cases. 

The dispersion relations plotted throughout this paper apply, strictly speaking, only 
to the on- but not off-shell improved actions obtained in the third step of our procedure. 
The reader might wonder what happens to the dispersion relations, when, as instructed 
in step 3, we undo the change of variable, leading to off-shell improved actions and fields. 
The off-shell improved propagator in eq. (3.6) implies the dispersion relation 

Tr{n{p)n*(p)) Tr{n{p)n*{p)) (f + m{pf) = . (C.IO) 
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The factors from undoing field transformation give rise to additional high-lying ghosts, 
without affecting the other branches of the dispersion relation. Considering, for 
simplicity, the SW case, where we can choose = U = 1 — "^{f — rric), we have 

\Tr{Q{p)Q*{p)) = \Tr{nip)n*{p)) = (l + '-^m^f + C-^fp' . (C.ll) 

For zero mass and momentum the extra ghosts are at energies aoE{0) = arcsinh(4/r). 
If these ghosts are deemed to be not high enough, one can push them up even further 
by modifying the field transformation operators to read 

(C.12) 

for suitably large integer n. Obviously, similar remarks apply to the D234 actions. 

So far we have considered the dispersion relation in the absence of a background 
gauge field, finding that the coefficients of the scaling errors of our D234 actions are 
quite small. One might wonder if small coefficients prevail in the presence of a gauge 
field. 

The short answer is yes, because, as for vanishing gauge field, these errors ultimately 
arise from the truncation of the expansion of the the continuum in terms of lattice 
derivatives, which is a series with a finite radius of convergence and rapidly decreasing 
coefficients. 

The longer answer goes as follows. When wc discrctizc the continuTim action Mq 
obtained by a change of variable in sect. 3.1, the leading error is of even order, a^, a*^, 
. . . , and comes from the truncation of the p term; the truncation error of always 
comes with an additional factor of uq. 

The leading error of the general dispersion relation of a D234 action (at least the 
D234c(r) one, see below) is therefore simply that of the naive improved fermion action 

+ rric- Remembering that Vc^ = — ^aj^D^ + 0{a^) we find the position-space 
dispersion relation 

ml = ^2^1 + l^-F - ^E<^M - ^E^a^AD^atDl] + 0{a') . (C.13) 

H M HI/ 

The first two terms on the rhs give the continuum dispersion relation, the rest the 0{a'^) 
errors. Note that the commutator terms vanish in the absence of a gauge field. As 
promised, the coefficients of the error terms are small. 

For the "tuned" D234 actions of sect. 3 there are additional errors, e.g. there is a 
— j^UqAq correction to the action for the D234i(|) case. Again, the coefficients of these 
errors are quite small. 

For comparison, here is the dispersion relation of the SW action 

ml = EL>2 ^ 1^.^ ^ iT^^Df, + ^Y.i^^'^iD^,alDl] + 0{a'). (C.14) 



n 



zLr, 



rric) 
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Note that using these expansions it is trivial to calculate the finite-mass corrections 

to various terms in the action. For example, for the case of a small, constant, magnetic 
background field one finds that the energy of a "D234c particle" at zero velocity is: 

= ml(l + ^atmt) - 1 a-F (l + ^ a^m^) + 0{a') . (C.15) 

The first term shows the correction corresponding to the leading coefficient in table 1. 
The second shows that the correction to the hyper-fine splitting, although not as small 
as for the mass, is still small; on the one percent level as long as aorric < 0.5. This should 
be compared with the SW case, where the analogous calculation gives 

E^ = ml(l- ^almfj - 1 a-F (l - ^ almfj + 0{a^) , (C.16) 

showing that aQtric must be smaller by almost a factor of 4 for the hyper-fine splitting 
correction to be similarly small. 

Appendix D Tadpole Calculus 

A large, apparently often dominant, part of the contributions of standard lattice 
perturbation to a generic quantity are unphysical, in the sense that they are due 
to tadpole diagrams, which do not exist in standard ways of performing continuum 
perturbation theory. Lepage and Mackenzie [8] proposed a simple method to improve 
lattice operators and actions, i.e. make them more continuum-like. 

At tree level, their prescription amounts to replacing each link field U^{x) in an 
operator by Un{x)/un, where the number < ti^ < 1 should be defined such that, 
roughly speaking, 1 — represents the "tadpole part" of the link field. At higher 
orders, one divides each link by the non-perturbatively measured u^, and multiplies by 
its perturbative expansion. In other words, one reorganizes perturbation theory in a 
manner that sums up most of the tadpole contributions to all orders, separately from 
the physical contributions. For isotropic actions all are the same and a simple, gauge 
invariant prescription for the is given hy uq = {^Tr Un)^^^, where Un is the plaquette 
operator in a theory with SU(A/") gauge fields. 

This tadpole improvement prescription is the one used in the literature up to now, 
and a version for anisotropic lattices was used for exploratory simulations [11, 23]. We 
are presently also exploring other prescriptions [18]. 

In the rest of this section we just assume that the have been chosen by some 
scheme; their precise numerical values are irrelevant. For our fermionic actions the 
operators to be tadpole-improved are V^i^, and various products and sums thereof. 
Let T{A) denote the tadpole-improved version of some operator A. Note that T is a 
linear operator, i.e. 

T{A + \B) = T{A) + XT{B) for any A gC . (D.l) 
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Clearly, 



V„ 



m,) = ^ + ^ - 1) • (D-2) 

A subtlety arises for sufficiently complicated products of operators. Namely, the 

improvement prescription requires one to expand out all products in terms of C/^(a;)'s 
and use U^{x)U-fj^{x + /x) = 1 for "backtracking" products of link fields. This implies 
that in general 

T{AB) ^ T{A) T{B) . (D.3) 

In practice it will sometimes be more efficient to expand out a product of operators 
in terms of link fields — then no subtlety arises. In other cases, however, it will be 
more convenient or efficient to apply one tadpole-improved operator after another. If we 
refer to the replacement of a product of Vyij's and A^^'s by the product of their tadpole- 
improved versions as "naive tadpole improvement" , the simplest strategy in such a case 
is to first apply naive tadpole improvement and then correct the error by using the 
following formulas: 

r(v^A^) = T(v^)r(A^) 
r(v^) = r(v^)2 + ^^ 

T{Al) = TiA.f - ^6,TiA,) + (D.4) 



where 

We also remind the reader of the identities (B.7) and (B.8). Easy corrollaries of eqs. (D.4) 
are 



Tif) = Tiff + C£S,)T{f) + (D-6) 
which are useful for the tadpole improvement of the field redefinition operators J7 and 

n. 
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